CHANGING THE BRANCHING MECHANISM OF A CONTINUOUS 
STATE BRANCHING PROCESS USING IMMIGRATION 
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Abstract. We consider an initial population whose size evolves according to a continu- 
ous state branching process. Then we add to this process an immigration (with the same 
branching mechanism as the initial population), in such a way that the immigration rate is 
proportional to the whole population size. We prove this continuous state branching process 
with immigration proportional to its own size is itself a continuous state branching process. 
By considering the immigration as the apparition of a new type, this construction is a nat- 
ural way to model neutral mutation. It also provides in some sense a dual construction 
of the particular pruning at nodes of continuous state branching process introduced by the 
authors in a previous paper. For a critical or sub-critical quadratic branching mechanism, it 
is possible to explicitly compute some quantities of interest. For example, we compute the 
Laplace transform of the size of the initial population conditionally on the non extinction 
of the whole population with immigration. We also derive the probability of simultaneous 
extinction of the initial population and the whole population with immigration. 



1. Introduction 

We consider an initial Eve-population of type whose size evolves as a continuous state 
branching process (CB), = {Y^,t > 0), with branching mechanism ip^ defined by 



:i) V (A) = a^A + /3A^ + / 7r{di) 

'(0,oo) 



where a" G M, /5 > and vr is a Radon measure on (0, oo) such that J^^ oo)(l '^^^) ^('^^) 

< oo. 

See [9] for a definition of CB as limit of Galton- Watson processes. We assume that this 
population undergoes some irreversible mutations with constant rate, giving birth either to 
one individual of type 1 (with rate a), or to infinitely many offsprings of type 1 (with rate and 
mutant offsprings size described by a measure v). This second population of type 1 evolves 
according to the same branching mechanism as the Eve-population (i.e. the mutations are 
neutral). The population of type 1 undergoes also some mutations and gives birth to a 
population of type 2 with the same rules, and so on. 

If we loose track of the genealogy, the new population of type 1 can be seen as an immigra- 
tion process with rate proportional to the size of the Eve-population, the population of type 
2 is an immigration process with rate proportional to the size of the population of type 1, 
and so on. We are interested in the law of the total population size X = {Xt,t > 0), which is 
a CB with immigration (CBI) rate proportional to its own size. If the mutations are neutral, 
we expect X to be a CB. This is indeed the case : if ip^ is the branching mechanism of the 
Eve-population and 

(/>(A) = a\+ [ u{dx) (l - e-^A 

JiO,+oo) ^ ' 
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is the immigration mechanism, then the total population size is a CB with branching mech- 
anism ijj = — (j)^ see Theorem 13.31 

Another approach is to associate with critical or sub-critical CBs a genealogical structure, 
i.e. an infinite continuous random tree (CRT), see [lOj or [5j. In that context, each individual 
of the CB can be followed during its lifetime and mutations can be added as marks on its 
lineage. Pruning the CRT associated with the total population (of branching mechanism 
ip = ip^ — (p'j at these marks allow to recover the Eve-population from the total population. 
This construction has been used in [2] and the construction given here via immigration 
proportional to the size of the population can be seen as the dual of the pruning construction 
of [2j , see Section S] and more precisely Corollary 14. 2[ However [2j considers only the case 
where the branching mechanism of X is given by a shift of the branching mechanism of the 
Eve-population. We shall give in a forthcoming paper [3] a more general pruning procedure 
which will correspond to the general proportional immigration presented here. 

Natural questions then arise from a population genetics point of view, where only the 
whole population Xt is observed at time t. In order to compute some quantities related 
to the Eve-population, given the total population, we compute the joint law of the Eve- 
population and the whole population at a given time: (1^^, X^). For quadratic critical or sub- 
critical branching mechanism, we provide explicit Laplace transform of the joint distribution 
of (1^*^, Xt). In particular, we compute f{Y^ = 0\Xt > 0), the probability for the Eve-type to 
have disappeared at time t, conditionally on the survival of the total population at time t, see 
Remark 15.31 We also compute the Laplace transform of conditionally on the population 
to never be extinct, see Proposition 15.61 In Lemma 15.51 we compute the probability of 
simultaneous extinction of the Eve-population and the whole population, in other words, 
the probability for the last individual alive to have undergone no mutation. The techniques 
used here didn't lead us to an explicit formula but for the quadratic branching mechanism. 
For the general critical or sub-critical case, we use, in [Tj, a Williams' decomposition of the 
genealogical tree to give a very simple formula for the probability of simultaneous extinction 
of the Eve-population and the whole population. 

In the particular case of CB with quadratic branching mechanism {ip{u) = fdu^, f3 > 0), 
similar results are given in [14J (using genealogical structure for CB) and in [15j (using 
a decomposition of Bessel bridges from [T2]). In the critical ('ip'{0~^) = 0) or sub-critical 
(V''(0"'") > 0) case one could have used the genealogical process associated to CB introduced 
by [To] and to CBI developed by [8] to prove the present result. This presentation would have 
been more natural in view of the pruning method used in [2]. Our choice not to rely on this 
presentation was motivated by the possibility to consider super-critical cases {ip'{0'^) < 0). 

The paper is organized as follows: In Section [21 we recall some well known facts on CB and 
CBI. In Section [31 we built a CBI X whose branching mechanism is and immigration rate 
at time t proportional to Xt and prove that this process is again a CB. We give in Section 
[3] some links with the pruning at nodes of CB introduced in [2j . Eventually, we compute 
the joint law of the Eve-population and the whole population in Section [5l as well as some 
related quantities. 



2. CB AND CB WITH IMMIGRATION 



The results from this section can be found in [7J (see also [TT] for a survey on CB and CBI, 
and the references therein). Let ^ he a branching mechanism of a CB: for A > 0, 



(0,oo) 



TT{d£) 



-1 + A^l 



{^<1} 



CSBP WITH IMMIGRATION 



3 



where a € M, /3 > and vr is a Radon measure on (0, oo) such that /^^ ^-,(1 A i'^) TT{di) < oo. 

Notice ip is smooth on (0, oo) and convex. We have ^''(O'^) G [— oo,+oo), and ip'{0^) = —oo 
if and only if Jj-^ £ 7r{di) = oo. In order to consider only conservative CB, we shall also 
assume that for all e > 

(3) / T—r^ du = oo. 



\^{n)\ 



Notice that ■(/''(O^) > — oo implies 



2.1. CB. Let be the law of a CB Z = {Zt,t > 0) started at x > and with branching 
mechanism ^/J. The process Z is a Feller process and thus cad-lag. Thanks to ([3]), the process 
is conservative, that is a.s. for all t > 0, Zt < +oo. For every A > 0, for every t > 0, we have 



(4) E,. 
where the function u is the unique non-negative solution of 

(5) u{t,X)+[ tp{u{s,\))ds = X, A>0, t>0. 



This equation is equivalent to 

dr 

(6) / -^=t A>0, t>0. 

Juit,\) i^ir) 

The process Z is infinitely divisible. Let Q be its canonical measure. The c-finite measure 
Q is defined on the set of cad-lag functions. Intuitively, it gives the "distribution" of the size 
process for a population generated by an infinitesimal individual. In particular, under P^, Z 
is distributed as Yli^i where X^jg/ ^z^ is a Poisson point measure with intensity xQ{dZ). 
Thus, for any non-negative measurable function F defined on the set of cad-lag functions, we 
have the following exponential formula for Poisson point measure 

E,,[e-S.e.^{^0] =exp(-xQ[l-e-^(^)]) . 



The CB is called critical (resp. super-critical, resp. sub-critical) if -(/"'(O^) = (resp. 
^^'(0+) < 0, resp. V'(0+) > 0). 

We shall need inhomogeneous notation. For t < 0, we set Zt = 0. Let Vx,t denote the law 
of [Zg-t, s G M) under P^;, and let Qt be the distribution of (Zg-t, s G M) under Q. 

For fi a positive measure on R, we set H^^ = sup{r G M;/i([r, oo)) > 0} the maximal 
element of its support. 

Proposition 2.1. Let fi be a finite positive measure on M with support bounded from above 
(i.e. is finite). Then we have for a// s G M, x > 0, 



(7) E, 



^-xw{s) 



where the function w is a measurable locally bounded non-negative solution of the equation 



(8) w{s)+ '4){w{r))dr = / /i((ir), s < and u;(s) = 0, s > . 

Js J[s,oo) 

lfip'{0^) > —oo or if fj,{{H'^''}) > 0, then ^ has a unique measurable locally bounded non- 
negative solution. 
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This result is well known for critical and sub-critical branching mechanism (see eg. [6]). 
As we didn't find a reference for super-critical branching mechanism, we give a short proof 
of this Proposition. 

Proof. Let n > 1. We set Z^'"'^'' = for t G [i/2", (i + l)/2''). Using that Z is cad-lag, 

we get a.s. lim zi"'^''^ = Zt-s for all t, s € M. Since the process Z is finite, we get by the 

?i— >oo 

dominated convergence theorem a.s. for all s € M 

Zr^s Kdr) = lim / zjP^'' fi{dr). 
Using the Markov property of Z, we get that 



Ex 



^ 1 



where it;^"^ is the unique non-negative solution of 

/■([/f'^2"] + l)/2" r 

\s)+ / ij{w^''\r)) dr = / fx{dr), 

Js J\k/2"-.oo) 



with k s.t. k/l"^ < s < {k + l)/2". 

Let T > Rf" + 1. Notice that for all s € [-T,T], we have / Z^""^'' n{dr) < sup{Zt,t € 
[0,2T]}^l{[-T,H^']) < oo a.s. Let C be defined by e"^ = E^.[e-™P^^*'*eI°'2^1>'^([-^-f^''l)]. 
Notice C < oo. This implies that for all n > 1, s G [—T,T], 

< w(")(s) < C < oo. 

By dominated convergence theorem, w^^\s) converges to w{s) = — log ^Ei[e"-^^''"*^^'^'''* 

which lies in [0, C], for all s E [— T, T]. By dominated convergence theorem, we deduce that w 
solves dSD. Since T is arbitrary, the Proposition is proved but for the uniqueness of solutions 

of dSD. 

If 'iIj'(0~^) > — oo, then tp is locally Lipschitz. This implies there exists a unique locally 
bounded non-negative solution of dH). 

If ^'(O"*") = — oo, and fi{{H'^}) > 0, we get that f Z^-s fJ-idr) > aZnt^^s, where a = 
fi{{H^}) > 0. This implies that w{s) > u{H^ — s,a) > for s G M. The function u{-,a) is 
strictly positive on M+ because of condition ([3]) and equation Since is locally Lipschitz 
on (0, oo), we deduce there exists a unique locally bounded non-negative solution of ([8|). □ 

2.2. CBI. Let x > 0, a > 0, be a Radon measure on (0, oo) such that J^^ ^^(1 Ax) i>{dx) < 
oo. Let denote the set of non- negative measurable functions defined on M. Let h E 
be locally bounded. We consider the following independent processes. 

• "^i^i Sp^x\z^j ^ Poisson point measure with intensity /i(t)l{^>o}(it i^{dx) Fx^t{dZ). 

• Z, distributed according to P^. 

• J2jeJ^p Z3'> ^ Poisson point measure with intensity Q/i(t)l{t>o}dt Qt{dZ). 

For t G M, let Yt = Zt + Y^iai Zt + Y^j^J ^ [0, oo]. We say Y = [Yt, t > 0) is a continuous 
state branching process with immigration (CBI) started at x, whose branching mechanism 
is ip and immigration is characterized with {h, cj)) where the immigration mechanism, <j), is 
defined by 

(9) (t>{X) = aX+ [ z^(dx)(l -e"^^), A > 0, 
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where a > 0, and is a Radon measure on (0, oo) such that J^q oo)^^ ^ "^^ i'{dx) < oo. 

One gets y is a conservative Hunt process when h is constant, see 0. Notice that y is a 
non-homogeneous Markov processes. We also have Yq = x, and Yj = for t < 0. 

Using Poisson point measure property, one can construct on the same probabihty space 
two CBI, Y^ and Y'^, with same branching process ■0, same starting point and immigration 
characterized by {h^,4>) and (/i^, </>) such that Y"/ < Y^ for all t < T as soon as h^{t) < h?{t) 
for all t <T. We can apply this with = h and h'^{t) = sup{/i(s); s G [0,T]} for t G M and 
some T > 0, and use that Y^ is conservative (see [7j) to get that Y^ has a locally bounded 
version over [0,T]. Since T is arbitrary, this implies that any CBI has a locally bounded 
version. We shall work with this version. 

The following Lemma is a direct consequence of the exponential formula for Poisson point 
measures (see eg. [13], chap. XII). 

Lemma 2.2. Let fi be a finite positive measure on M with support bounded from above (i.e. 
is finite). We have for s G M.' 

p~^'^(*)~/o°° h{t)<p(w(s+t))dt 

where the function w is defined by 



(10) 



E 



3. State dependent immigration 

3.1. Induction formula. Let {xk,k G N) be a sequence of non-negative real numbers. Let 
Y^ be a CB with branching mechanism tp^, defined by ([T]), starting at xq. We shall assume 
that Y^ is conservative, that is condition 1^ holds for tp^. We construct by induction y", 
n > 1, as the CBI started at x„, with branching mechanism ^jj^ and immigration characterized 
by (y"~\ (/>), with (f) given by ([9]). 

Lemma 3.1. Let (^fc. A; G N) be a family of finite measures on R with support bounded from 
above. We have for a// n G N, s G M, 

p-Efe = 0^n-fe"'fe"^(«) 



E 



where Wq"^ is defined by ^ with n replaced by ^n, and for k>l, is defined by ^ with 
fi replaced by Hn-k{dr) + (l){wk-i{r)) dr. In particular, Wk is a locally bounded non-negative 
solution of the equation 

11) wis)+ I ^\w{r))dr=l fin^kidr)+ j Hwt'-iir)) dr, sG 

(Notice Wk{s) = for s> max{iJ''fc' , /c' G {0, . . . , k}}.) 



y°,...,y"-i 



Proof. This is a consequence of the computation of E 

using Proposition 12. 1[ This also implies that pTj) holds. Then, by induction, one deduces 
from (jlip that Wk is locally bounded. □ 



3.2. Convergence of the total mass process. We consider the sequence (Y"-,n > 0) 
defined in the previous Section with xq = x > and 3;„ = for n > 1. We set = Ylk=o 
for t G M. Let Xt be the increasing limit of X" as n ^ -|-oo, for all t G M. We have 
Xt G [0, +CX3]. We call X = {Xt,t G M) a CBI with branching mechanism ip^ and immigration 
process {X, (p). We set ip = ip^ — (p. 
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Remark 3.2. For A > 0, we have 

V'°(A) - (piX) = (a^-a- [ i i^{di)) X + pX^+ [ {t^{M) + v{dl)) [e"^^ -1 + A£lo<n 

^ J(0,1] ^ i(0,oo) L 

This gives that iIj = ip^ — (f) \s a. branching mechanism. 

The process describes the size process of the Eve-population, the size process of 
the population of mutants born from the Eve-population , Y'^ the size process of the 
population of mutants born from mutant population Y^, and so on. The size process of the 
total population is given hy X = "^k^oY^. In neutral mutation case, it is natural to assume 

that all the processes Y^ have the same branching mechanism. Since we assume Xk = for 
all A; > 1, this means only the Eve-population is present at time 0. 

Theorem 3.3. We assume that ip is conservative, i.e. satisfies The process X, which is 
a CBI with branching mechanism iJj^ and immigration process (X, (j)), is a CB with branching 
mechanism ip = tp^ — (p. 

Remark 3.4. As a consequence of Theorem l3.3|, X is a Markov process. Notice that {Y^, . . . , Y"' 
is also Markov but not {X^,t > 0) for n > 1. 

Proof. Let be a finite measure on M with support bounded from above (i.e. < oo). We 
shall assume that ^{{H^}) = a > 0. 

We keep the notations of Lemma 13.11 with fj.^ = fj,. In particular we see from (jlip that 

(n) 

w) does not depend on n. We shall denote it by Wk- By monotone convergence, we have 



E 



hm E 

n— >oo 



lim e-^"'"^"). 



where the limits are non-increasing. This implies that {wn, n > 0) increases to a non- negative 
function w^o- By monotone convergence theorem (for V'°(^(^))l{«)n{r)>o} ^'^'^ the 
integral with (p) and dominated convergence theorem (for V'°(^(''))l{ti;„{r)<o} dr)^ we 
deduce from pT]) . that w^o solves w{s) = d iov s > H^^ and 



0/ 



(12) w{s) + / i^^{w{r))dr = / ^(dr) + / (p{w{r)) dr, s < if^. 

Js J[s,oo) Js 

Notice that Woo{s) € [0, oo] and the two sides of the previous equality may be infinite. 

Thanks to Proposition 12. H and since ip^ — cp is a branching mechanism (see Remark | 
there exists a unique locally bounded non- negative solution of (jl2p . which we shall call w. 
Therefore to prove that Woo = w, it is enough to check that Woo is locally bounded. This will 
be the case if we check that Woo < w. In particular, we get Woo = w, if we can prove that 
Wn ^ w for all n S N. We shall prove this by induction. 

We consider the measure ii^{dr) = ij.{dr) + (p{w{r))l^r<Hi^} dr. Notice Z/^" = and 



^^{{H^""}) = /u({F^}) = a > 0. We define ivq by 



e 



The function wq is a locally bounded non-negative function which solves 



0/ 



w{s) + I ir{w{r))dr = / ^(dr) -|- / (p{w{r)) dr, s < H^. 

:,00) 
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Thanks to Proposition l2.lt wq is unique. Since w solves the same equation, we deduce that 
ivQ = w. We also have 



^-xwo{s) _ g 



> E 



This implies that wq < wq = w. 

Assume we proved that Wn-i < w for some n > 1. Then we can consider the measure 
^"(dr) = ^l{dr) + [(t){w{r)) - (/.(u;„_i(r))]l{,.<HM} dr. Notice H^'" = H'' and = 
a > 0. Recall x = xq > and Xk = ior k > 1. We define Wn by 



-xw„{s) _ Jg 



The function Wn is a locally bounded non-negative function which solves for s < 



wis) + 



ip^{w{r))dr 



s,oo) 



fi^idr) + 



(m„_i(r)) dr 



,oo) 



ld{dr) + 



(j){w{r)) dr. 



Thanks to Proposition 12. H Wn is unique. Since w solves the same equation, we deduce 
that Wr, = w. We also have 



> E 



This implies that Wn < w. Therefore, this holds for all n > 0, which according to our previous 
remark entails that Woo = w. 

By taking ^i{dr) = Y.k=i ^^kk^^) for G N*, Ai, . . . , e [0, oo) and < ti < . . . < 
tx, we deduce that X has the same finite marginals distribution as a CB with branching 
mechanism if:^ — (p. Hence X is a CB with branching mechanism — (p. 

□ 



4. The dual to the pruning at node 

For S M, we consider the group of operators {Tg,6 E M) on the set of real measurable 
functions defined by 

Te{f){-) = f{9 + .)-f{e). 

Let ijj^ be given by with Levy measure vr. Using the previous Section, for > 0, we 
can give a probabilistic interpretation to T^Q^ip^) as a branching mechanism of a CBI with 
proportional immigration. Let Oq = sup{0 > 0; /^-^ eP^ Tr{d£) < oo}. Notice that = 

if (0+) = — oo, as (0+) = —oo is equivalent to f^-^ i ■K{di) = +co. We assume 0o > 
and we set G = (0,6*0] if /(i 00) ^('^^) < 00 and 6 = (0,6lo) otherwise. Let 6* G 6. We 
define 

cpeW = 2peX + / (e^^' -1)(1 - e"^^) Tr{dx). 

It is straightforward to check that T_q{iIj^) = — (pQ and that (pg is an immigration mecha- 
nism. Notice that for 6 < 6q, we have T_g(Tp^y{0^) > — 00, that is T_g{ip^) is a conservative 
branching mechanism. 

The next Corollary is a direct consequence of the previous Section. 
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Corollary 4.1. Let 9 £ Q. If 9 = 9q assume furthermore that T^0g{'ip^) is conservative. A 
CB I process X with branching mechanism ip and immigration {X,(/)g) is a CB with branching 
mechanism T_g{'ip). 

On the other end, for 9 > 0, Tg{ip) can be seen as the branching mechanism of a pruned 
CB. The following informal presentation relies on the pruning procedure developed in [2j. 
Let us consider a CRT associated with a critical or sub-critical branching mechanism ip with 
no Brownian part, which we shall write in the following form: 

^(A) = QiA+/ (e-^^ - 1 - Ar)7r((ir), 

i(0,+oo) 

with J^-^ r Tr{dr) < oo, ai = a + J^-^ r Ti{dr) > and J^q rTr{dr) = +oo. In that case, 
the CB process with branching mechanism has no diffusion part (/? = 0) and increases only 
by positive jumps. Let us recall that a CRT can be coded by the so called height process 
H = {Ht, t > 0), see [5]. Intuitively, for the individual t > 0, Ht > represents its generation. 
The individual t is called an ancestor of s if Ht = mm{Hu,u € [s A t,s V t]}, and we shall 
write s )^ t. Informally, for t fixed, the "size" of the population at generation a > of 
all individuals r < t is given by the local time of H at level a up to time t, say. For 
the CRT associated with the branching mechanism ip, the process L = {L^^,a > 0), where 
Tx = inf{t > 0, Lj > x} is a CB with branching mechanism ^. The height process codes for 
the genealogy of the CB process L. 

An individual t is called a node of the CRT if the height process corresponding to its 
descendants, {Hg — Ht, s )p t) has a positive local time at level 0, say A^. (If t < Tx, then At 
corresponds to a jump of the CB process L at level Ht; reciprocally to a jump A of the CB 
process L at level a there corresponds an individual t < Tx such that Ht = a and A^ = A.) 
Intuitively At corresponds to the "size" of the offspring population of individual t. Let 9 > 
be fixed. A node t of size A^ is marked with probability 1 — e~^'^' , independently of the other 
nodes. To prune the CRT, we just remove all individuals who have a marked ancestor. The 
height process of the pruned CRT is then given by H^ = {Hc^ ,t > 0) where C is the inverse 
of the Lebesgue measure of the set of individuals whose ancestors have no mark: 

Ct = inf{ro > 0; / l{\/s,r)ps,s is not marked} dr > t} . 

Jo 

Theorem 6.1 in [2] shows that this pruned CRT is itself a CRT associated with the branching 
mechanism Tq[iP). 

By looking at the local time of the pruned process, we get a nice construction of a CB 
process of branching mechanism Tg^tp), which we shall called a pruned CB with intensity 
^ > 0, from a CB process of branching mechanism tp. Notice this construction was done 
under the assumption that (3 = (see also [4| when /3 > and vr = 0). The general pruning 
procedure in the general case /3 > and vr 7^ will be presented in a forthcoming paper [3]. 

In a certain sense the immigration is the dual to the pruning at node: to build a CB 
process of branching mechanism ip from a CB process of branching mechanism Tg^ip), with 
9 > 0, one has to add an immigration at time t which rate is proportional to the size of the 
population at time t and immigration mechanism (pg defined by: 

(^g(A) = Tg{iP){\) - ^(A) = 2I39\ + [ (1 - e-^^')(l - e"^^) Tr{dx), for A > 0. 

J(0,oo) 

In other words, we get the following result, whose first part comes from Theorem 6.1 in p]. 
As in |2], we assume only for the next Corollary that /? = and i ^{dtj = +00. 
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Corollary 4.2. Let X be a critical or sub-critical CB process with branching mechanism ip. 
Let X(^) be the pruned CB of X with intensity 9 > : X^^^ is a CB process with branching 
mechanism To^ip). The CBI process, X, with branching mechanism TQ{ij)) and immigration 
{X,(pg) is distributed as X. 



5. Application : law of the initial process 



We consider a population whose size evolves as X = {Xt,t > 0), a CB with branching 
mechanism ip given by ([2]). We assume ip satisfies the hypothesis of Section[2l This population 
undergoes some irreversible mutations with constant rate. Each mutation produces a new 
type of individuals. In the critical or sub-critical quadratic case (vr = 0) this corresponds to 
the limit of the Wright-Fisher model, but for the fact that the "size" of the population is not 
constant but is distributed 3jS Si CB. 

We assume the population at time has the same original Eve-type. We are interested in 
the law of = {Y^'^t > 0), the "size" of the sub-population with the original type knowing 
the size of the whole population. In particular, we shall compute f{Y^' = 0\Xt > 0), the 
probability for the Eve-type to have disappeared, conditionally on the survival of the total 
population at time t. 

We shall assume Y^ is a CB with branching mechanism and X is the CBI with immi- 
gration (X, (/>), with = -0^ — considered in Section [3.21 Thus, we model the mutations 
by an immigration process with rate proportional to the size of the population. 

The joint law of {Xt, Y^) can be easily characterized by the following Lemma. 



Lemma 5.1. Let t > 0, Ai, A2 G 



We assume Xq = Yr^ 



E 







x > 0. We have 



where {w,w*) is the unique measurable non-negative solution on (— oo,t] of 
w{s) -\- I ip'^[w{r))dr = \i + \2+ [ (l){w* {r))dr, 



Proof. Recall notation of Section 
apply Lemma l3. II with 



w*{s)-\- I tp^w* {r)^dr = Xi. 

In particular xq = x and x„ = for all n > 1. Let us 



fio{dr) = (Ai + X2)6t{dr), 
fj-kidr) = Xi5t{dr) for A; > 1. 



We get 



where for s <t, 



E 



ij\w^^\r))dT = Ai, 
V'°(u'["^(r))(ir = Ai 



(p(^w^}^{r)^dr for 1 < k < n — 1, 



I V°(t«l"^(r))(ir = A1 + A2+ / </'(«;2i(r))dr. 
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We let n goes to infinity and use similar arguments as in the proof of Theorem 13.31 to get the 
result. □ 

Some more explicit computations can be made in the case of quadratic branching mecha- 
nism (see also p[5] when a = 0). Let a > and 9 > and set 

tpiu) =au + u^, -0°^ = T0{ij){u) = {a + 2e)u + . 

The CB which models the total population is critical (q = 0) or sub-critical (a > 0). The 
immigration mechanism is (/'(u) = 7p^(u) — ip{u) = 26u. 
We set b= {a + 29) and for t > 0, 



(13) h{t) = <j 

Proposition 5.2. Let t > 0, Ai, A2 G 



l + Ai^^l^ ifa>0, 
1 + Ait ifa = 0. 

M+. We have 



E 



-XiXt-X2Yf 



-xvo{t) 



where 



voit) = e-'' h{t) 



f e-^" /i(r)-2dr ) + Ai e""* /i(t)-^ 
A2 Jo J 



Proof. By the previous lemma, we have 



E 



-xw{0) 



(14) 

where for s < t, 

(15) w{s) + I w{r){w{r) +b)dr = Xi + X2 + 2e f w*{r)dr, 

J s J s 

w*{s)+ / w* {r)[w* (r) + a)dr = Xi. 

J s 

The last equation is equivalent to 

(16) {w*)' -w*{w* +a) = ^ on(-oo,t], w*{t)=Xi. 
The function z* := — is thus the unique solution of 



{z*)' + az* + I = Q on(-oo,t], z*{t) 



Ai' 



If a > 0, this leads to 



z*{s) = -\ e 
a 

1 



a{t-s) 



a{t-s) 



If a = 0, we have z*{s) = t — s + —. We get 

Ai 

(17) w*{s) = h'{t - s)h{t - sy^ = Ai e-"(*-^) h{t - s)-\ 

Equation (fT5l) is equivalent to 

w' — w{w + b) = —26w* on (—00, t] , w^t) = Ai -|- A2. 
Set y = w — w* and use the differential equation (fTB|) . to get that y solves 
y' -y^ -y{2w* + b)=0 on(-oo,t], y{t) = X2. 
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Then the function z := 1/y is the unique solution of 

z' + {2w* +b)z + l = on(-oo,t], z{t) 
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1 

a; 



One solution of the homogeneous differential equation z'q = —{2w*+b)zo is zq{s) = e^^^~'^^ h{t- 
s)^. Looking for solutions of the form z{s) = C{s)zq{s) gives 



z{s) = zq{s) + zq{u) ^du^ . 



We conclude using (|M|) and w = w* + z . 



□ 



Remark 5.3. We can compute the conditional probability of the non extinction of the Eve- 
population: ^{Y^ > 0\Xt > 0). However, this computation can be done without the joint 
law of {Xt,Y^°) as 



F{YI^ >0\Xt>0) 



F{Yt° >0,Xt>0) _ F{Yt^ > 0) _ 1 - F{Yl> = 0) 
F{Xt > 0) ~ F{Xt > 0) ~ 1 - F(Xt = 0) '■ 



with F{Xt = 0) = lim E[e"^i^*] 

Ai— >oo 



(18) 



g{a,t) 



^-xg{a,t) 1 p^y^o = Q) = lim E[e-^2F0] 

A2— +00 



^ ifa>0, 
t if a = 0. 



The same kind of computation allows also to compute the joint law at different times. 
Proposition 5.4. Let < u < t, Ai, A2 € M+. We have 



E 



-xv\ {u,t) 



where 

vi{u,t) =e-^*h{ty^ 



A^ 



+ / e-^*- /i(r)-2dr + Ai e""* 



t—u 



Proof. Recall notation of Section [3.2[ In particular xq = x and = for all n > 1. Let us 
apply Lemma |3. II with 

/Uo = Ai(5t + A25m, 
/ijt = Ai(^t for /c > 1. 
Let n goes to infinity as in the proof of Lemma 15.11 to get that 



(19) 



E 



where {w,w*) is the unique non-negative solution on (— oo,t] of 



(20) 



{s)+ ij^{w{r))dr = \i + \2l{s<u}+ (p{w*{r))dr. 



w 



*(s)+ / ij{w*{r))dr = \i 



Notice w* is still given by (jl7p . For s > u, we have w{s) = w*{s) and, for s < u, Equation 
(pOj) is equivalent to 

w' — w{w + h) = —29w* on(— oo,m], w{t) = w* {u) + \2- 
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From the proof of Proposition 15.21 we get 



■w{s) — w*{s) 
We conclude using ([19 



h{t - s) 



-b{t~u) h{t-uy 



+ 



□ 



At this stage, we can give the joint distribution of the extinction time of X, tx = inf {t > 
0-Xt = 0}, and of y°, ryo = inf{t > 0;Yt° = 0}. For -u < t, we have P(rx < t,ryo < u) = 
hmAi^oo, Aa^oo exp -xvi{u, t) that is 



{tx < ^1 Tyo < u) = exp —X ( e 



J t—u 



We can compute the probabihty of simultaneous extinction of the Eve-population and the 
whole population, see also proposition 5 in [15], where a = 0. In [Ij, using different techniques 
we derive this formula for general critical or sub-critical branching mechanisms. 



Lemma 5.5. We have P(ryo = tx\tx = t) = e 



-2et 



Proof. We have Firyo = tx\tx = t) = 1 ■ — - — ^ 



-2et 



□ 



We can deduce from the latter Proposition the law of 1^ conditionally on the non-extinction 
of the whole population. We set 

A{b,u) = ^e'^+g{h,u). 
Proposition 5.6. Let u > 0, A2 G M+. We have 



lim E 



where 



G{a,u) = — e^'^ g{a,u) + 
A2 



^xA{h,u)-^ [l- A{h,u)-^G{a,u)) , 
2dig{b, u) if a = 0. 



Proof. We have 



E 



Xt>0 



E 



-X2YS 



E 



{Xt=0} 



nxt > 0) 



Using Proposition [5 

E 



u 1 



{Xt=0} 



= Mm E 













with vi{u,t) = \imx-^^+ooVi{u,t). 
Definition (jlSp implies 



X2g{a,t - m)2 



vi{u,t) 



+ e 



bt 



t-u 9{a,ry 



+ e-"*5(a,t) 



-1 



Performing an asymptotic expansion of vi as t goes to 00 leads to the result. 



□ 
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